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We present a precision calculation of the transverse-momentum and invariant-mass distributions 
for supersymmetric particle pair production at hadron colliders, focusing on Drell-Yan like slepton 
pair and slepton-sneutrino associated production at the CERN Large Hadron Collider. We imple- 
' ment the joint resummation formalism at the next-to-leading logarithmic accuracy with a process- 

' independent Sudakov form factor, thus ensuring a universal description of soft-gluon emission, and 

, consistently match the obtained result with the pure perturbative result at the first order in the 

^Sj ' strong coupling constant, i.e. at 0[as). We also implement three different recent parameterizations 

_ ' of non-perturbative effects. Numerically, we give predictions for crc'ji production and compare the 

resummed cross section with the perturbative result. The dependence on unphysical scales is found 
. to be reduced, and non-perturbative contributions remain small. 

Q^ ' PACS numbers: 12.60. Jv,13.85.Ni,14.80.Ly 

^ ! 

^\ I. INTRODUCTION 

^: 

. One of the main tasks in the experimental programme of the CERN Large Hadron Collider (LHC) is to perform an 
extensive and conclusive search of the supersymmetric (SUSY) partners of the Standard Model (SM) particles pre- 
\ dieted by the Minimal Supersymmetric Standard Model [J, [1] . Scalar leptons are among the lightest supersymmetric 
^ ■ particles in many SUSY- breaking scenarios Presently, the experimental (lower) limits on electron, muon, and 

\ tau slepton masses are 73 GeV, 94 GeV, and 81.9 GeV, respectively Since sleptons often decay into the corre- 
: sponding SM partner and the lightest stable SUSY particle, the distinctive signature at hadron colliders will consist 
in a highly energetic lepton pair and associated missing energy. 

The leading order (LO) cross section for the production of non- mixing slepton pairs has been calculated in @, 0,11,0], 
On while the mixing between the interaction eigenstates was included in [HJ]. The next-to-leading order (NLO) QCD 
, corrections have been calculated in iJJj], and the full SUSY-QCD corrections with non-mixing squarks in the loops 
' have been added in [T^ . Recently, an accurate calculation of the transverse-momentum [qx) spectrum including 
soft-gluon resummation at the next-to-leading logarithmic (NLL) accuracy has been performed [13|, allowing for the 
reconstruction of the mass and the determination of the spin of the produced particles by means of the Cambridge 
(s)transverse mass variable [H, 111 and for dist inguishing thus the SUSY signal from the SM background, mainly 



o 



X 



due to WW and tt production 



161 . |17| • Very recently, the mixing effects relevant for the squarks appearing in the 
loops have been investigated at NLO, and the threshold-enhanced contributions have been computed at NLL p^ . 
The numerical results show a stabilization of the perturbative results through a considerable reduction of the scale 
dependence and a modest increase with respect to the NLO cross section. 

Since the dynamical origin of the enhanced contributions is the same both in transverse-momentum and threshold 
resummations, i.e. the soft-gluon emission by the initial state, it would be desirable to have a formalism capable to 
handle at the same time the soft-gluon contributions in both the delicate kinematical regions, qx ^ M and ^ s, 
M being the slepton pair invariant-mass and s the partonic centre-of-mass energy. This joint resummation formalism 
has been developed in the last eight years [l^, [2^ . The exponentiation of the singular terms in the Mellin (N) and 
impact-parameter (6) space has been proven, and a consistent method to perform the inverse transforms, avoiding the 
Landau pol e and the singularities of the parton distribution functions, has been introduced. Applications to prompt- 
photon [2l[ , electroweak boson [22| , Higgs boson [IJ , and heavy-quark pair f23| production at hadron colliders have 
exhibited substantial effects of joint resummation on the differential cross sections. 
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In this paper we apply the joint resummation formalism at the NLL level to the hadroproduction of slepton pairs 
at the LHC, thus completing our programme (started in Ref. [l^ and continued in Ref. of providing the first 
precision calculations including soft-gluon resummation for slepton pair production at hadron colliders. In Sec. [ill 
we briefly review the theoretical formalism of joint rcsmnmation following Refs. (20l. [2^. We reorganize the terms of 
the resummed formula in a similar way as it was done for transverse-momentum resummation in ^25.] . The inverse 
transforms from the Mellin and impact-parameter spaces and the matching of the resummed result with the fixed- 
order perturbative results are discussed in Sec. IIIII Sec. IIVI is devoted to phenomenological predictions for the LHC, 
together with a comparison of the three types of resummation (transverse-momentum, threshold, and joint), showing 
their impact on the gy-spectrum and on the invariant-mass distribution. Our results are summarized in Sec. |Vl 



II. JOINT RESUMMATION AT THE NEXT-TO-LEADING LOGARITHMIC ORDER 



We consider the hard scattering process 

haiPa) hbipb) ^ F{M, qr) + X, (1) 

where is a generic system of colourless particles, such as a Higgs boson or a Drell-Yan (s)lepton pair, AI is the 
invariant mass of the final state F, and qx is its transverse momentum. Thanks to the QCD factorization theorem, 
the unpolarized hadronic cross section 

dV V- . ^ . . ^ d^aafc 



dAP Aqf, 

^ a,b 



can be written as the convolution of the relevant partonic cross section cTab with the universal distribution functions 
fa,b/ha b '-'^ partons a, 6 inside the hadrons /la.fc, which depend on the longitudinal momentum fractions of the two 
partons Xa,b and on the unphysical factorization scale fip- The partonic scattering cross section depends on the strong 
coupling constant a^, the unphysical renormalization and factorization scales and hf, and on the scaling variable 
z = M'^/s, where s = XaXbS and S = (pa + Pb)'^ are the partonic and hadronic centre-of-mass energies, respectively. 
The lower limits for the integration over the longitudinal momentum fractions contain the quantity r = M'^/S, which 
approaches the value t — 1 when the process is close to the hadronic threshold AP ~ S. In Mellin iV-space, the 
hadronic cross section naturally factorizes 

E£ ^^r-^ fa/kAN+l;iiF)fb/hAN + l;f^F)^f^iN;as,fiR,^iF), (3) 



dAP dq 



2 



,.b 



where the contour C in the complex TV-space will be specified in Sec. IIIII and the A^-moments of the various quantities 
are defined according to the Mellin transform 



1 

N-l 



F{N) ^ / dxx'^-'F{x) (4) 
Jo 

for X = Xa,b, z, T and F = fn/h^ b/hh^^^\'^\ respectively. The jointly resummed hadronic cross section in A'^-space can 
be written at NLL accuracy as 20l.l2^.[23t 

^2^(rcs) /■ (j2|-j 

dA'Pdq^ (^;«s,MJ?7/^f) = ^orfz Hcc{as,^lR) J —e'-^''^'^ Cc/ha{N,b; as, fiR, fip) 

X exp^El^'^'^{N,b;as,fiR) Cg/,,^, (A^, 6; a^, ^i?, ^f). (5) 

The indices c and c refer to the initial state of the lowest-order cross section a^g^ and can then only be qq or gg, since 
the final state F is assumed to be colourless. 

For slepton pair and slepton-sneutrino associated production at hadron colliders, 

haiPa) hb{pb) ^ k{pi)lf*{p2) + X, (6) 
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we have = {pi + P2) , = {v\t - Vrr) ■, and 



-(0) 



9M2 



^2 ^2^,, 



4a;H'(l-xn/)(l-m|/M2) 
2 



I il i ^ ^ lil j Z 



99' 



9M2 



(7) 



(8) 



where i, j denote slepton/sneutrino mass eigenstates with masses m^j , and mw are the masses of the electroweak 
gauge bosons, a is the electromagnetic fine structure constant, xw = sin^ 9w is the squared sine of the electroweak 
mixing angle, and the velocity /? is defined as 

(3 = ^1 + to4/M4 + m4/M4 - 2(to2/m2 + m^^/M^ + m| m^/M"). (9) 
The coupling strengths of the left- and right-handed (s)fermions to the electroweak vector bosons are given by 

{Lff'z,Rff'z} = {2Tf-2efXw)x6ff,, 
{^f^f'Z'^fJ'.z} {^ff'z Sl^*,Rff'z *}, 



{Lqq'W, Rqq'w} — {^2 Cw Vqq' , 0} , 
{Lq^q'Wi Rq.q'w} — {^qq'W S'^i S'^j^, 0} 



(10) 



where the weak isospin quantum numbers are Tj — ±1/2 for left-handed and = for right-handed (s)fermions, 

cw is the cosine of the electroweak mixing angle, and Vff are the CKM- matrix elements. The unitary matrices 5^ 
diagonalize the sfermion mass matrices, since in general the sfermion interaction eigenstates are not identical to the 
sfermion mass eigenstates (see App. |A|). 

The function Hce in Eq. ^ contains the hard virtual contributions and can be expanded perturbatively in powers 
of as, 



^.c>.A^.) = 1 + E(^) 



(11) 



The coefficients 

) = ^C,/b(A^;a.(M/x))C/6/a(A^;M/x,Mf)/a/hJ^ + l;A'F) (12) 

a.h 

and Cc/hb-i defined analogously, allow to evolve the parton distribution functions fa.b/ha b fro™ the unphysical factor- 
ization scale /Xi? to the physical scale M /x with the help of the QCD evolution operator 



^ dg2 

—lb/a{N]as{q)) 



(13) 



Ub/aiN;n,fJ,o) = exp 
and to include, at this scale, the fixed-order contributions 

oo 
n=l 

that become singular when qt ^ (but not when z ^ 1). The QCD evolution operator fulfils the differential equation 

dUb/a{N;fI,^lo) 

-^^[^ = 2^ Ub/c{N; pi, pio) lc/a{N; a,(/i)), (15) 
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where the anomalous dimensions ^c/a{N]as) are the iV- moments of the Altarelh-Parisi sphtting functions. The 
function 



N 



l + rib/N 



with b = bMe'"'/2 and iV ee A^e^^ 



(16) 



organizes the logarithms of b and N in joint resummation. Its exact form is constrained by the requirement that 
the leading and next-to- leading logarithms in b and N are correctly reproduced in the limits 6 ^ oo and N ^ oo, 
respectively. The choice of Eq. (jl6p with 77 = 1/4 avoids the introduction of sizeable subleading terms into perturbative 
expansions of the resummed cross section at a given order in as , which are not present in fixed-order calculations [2^ . 
The perturbative (PT) eikonal exponent 



(17) 



allows to resum soft radiation in the A-term, while the i?-term accounts for the difference between the eikonal approx- 
imation and the full partonic cross section in the threshold region, i.e. the flavour-conserving collinear contributions. 
In the large- limit, these coefficients are directly connected to the leading terms in the one-loop diagonal anomalous 
dimension calculated in the MS factorization scheme |2a| 



7c/c(A^; as) = -A^ias) InN - + 0{1/N). 

They can thus also be expressed as perturbative series in as, 

00 00 
A,(a.) = ^(^)"4") and Bc(as) = ^ (f)" ^' 

n— 1 n— 1 

Performing the integration in Eq. (jl7p . we obtain the form factor up to NLL, 



(») 



(18) 



(19) 



(20) 



with 



Ai^' 2A + ln(l-2A) 
To 

A^ 

f3'o 



X 

1 , 2 A In (1-2 A) 

2^-(^-^^)+ 1-2A 



/3o f^R 



A 



(2) 



2A 



1 - 2A 



+ ln(l - 2 A) 



B. 



(1) 



/3o 



ln(l - 2 A) 



and A = /Jq/tt as (/ii?) In x. The first two coefficients of the QCD /3-fmiction are 



^{UCA-'iTRNf) and (3^ ^ ^{17 C^- IOTrCaNj - GCpTRNf), 



(21) 



(22) 



Nf being the number of effectively massless quark flavours and Cp ~ 4/3, Ca ~ 3, and Tr = 1/2 the usual QCD 
colour factors. 

In order to explicitly factorize the dependence on the parameter x, it is possible to reorganize the resummation of 
the logarithms in analogy to the case of transverse-momentum resummation psi . [27j . The hadronic resummed cross 
section can then be written as 



dAf2 dgf. 



a.b 



dN 



y Hab^cc(^N;as, hr, ^J-F^ exp [^/^(Inx; ^s, M-r)] • 



(23) 
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The function Tiab^cs does not depend on the parameter x and contains all the terms that are constant in the limits 
6 ^ 00 or iV ^ oo, 



L n=l \ / 



(24) 



At 0{as), the coefficient H^^b-^cc given by 

^^aLc-c{N; I^f) - 6,aS,b Hil\f,n) + 6,a C^}l{N) + 5,b cl'/liN) + {s,ali}l{N) + 6,bl^^/^{N)) In ^. (25) 



The x-dependence appearing in the C-coefficient and in the evolution operator U of Eq. is factorized into the 
exponent Qc, which has the same form as -Ec^"""^ defined in Eq. (|17p except for the substitution 



B,{a,) ^ B,{N- as) = B,(a,) + 2/3(a,) ^''^ + 2 7,/,(iV; a,). 

d in Qfs 



(26) 



At NLL accuracy, Eq. (PO)) remains almost unchanged, since only the coefficient g^c^ of Eq. (pij) has to be slightly 
modified by 



B(i)_B(i)(Ar)^B(i)+27W(iV). 



(27) 



Although the first-order coefficients C^^^^{N) and H^^{^p) are in principle resummation-scheme dependent [27], this 
dependence cancels in the pcrturbative expression of Tiab^cc [25]. In the numerical code we developed for slepton 
pair production, we implement the Drell-Yan resummation scheme and take Hqq{as, ^jlr) = 1. The C-coefficients are 
then given by 



and 



' 2 (TV + 1) (TV + 2) 



1 



(28) 



III. INVERSE TRANSFORM AND MATCHING WITH THE PERTURBATIVE RESULT 

Once resummation has been achieved in N- and 6-space, inverse transforms have to be performed in order to get 
back to the physical spaces. Special attention has to be paid to the singularities in the resummed exponent, related 
to the divergent behaviour near x — exp[7r/(2/?oQ;s)], i.e. the Landau pole of the running strong coupling, and near 
b = —2N and b — —AN, where x = and infinity, respectively. The integration contours of the inverse transforms in 
the Mellin and impact parameter spaces must therefore avoid hitting any of these poles. 

The &— integration is performed by deforming the integration contour with a diversion into the complex 6-space 
PH . defining two integration branches 

b — {cos (fi ± i sin (p)t with < t < oo, (29) 

valid under the condition that the integrand decreases sufficiently rapidly for large values of \b\. The Bessel function 
Jo is replaced by two auxiliary functions hi^2iz,v) related to the Hankel functions 



hi(z,v) = / d6'e 



iz sin 6 



TT 



— tVTT 
iv-TT 



h2iz,v) = -3 / dee-" (30) 



Their sum is always hi{z, v) + h2{z, v) — 2 Jq{z), but they distinguish positive and negative phases of the 6-contour, 
being then associated with only one of the two branches defined in Eq. ((29)) . 

The inverse Mellin transform is performed following a contour inspired by the Minimal Prescription [2^ and the 
Principal Value Resummation (29j . where one again defines two branches 

N = C + ze'^"'' with < z < oo, 7r>(/)>^. (31) 
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The parameter C is chosen in such a way that all the singularities related to the iV-moments of the parton densities 
are to the left of the integration contour. It has to lie within the range < C < exp[7r/(2/?oas) — 1e\ in order to 
obtain convergent inverse transform integrals for any choice of (j) and if. 

A matching procedure of the NLL resummed cross section to the NLO result has to be performed in order to keep the 
full information contained in the fixed-order calculation and to avoid possible double-counting of the logarithmically 
enhanced contributions. A correct matching is achieved through the formula 



dM2 dg2 



dV^ 
dM2 dql 



Cn 



2ni 



bdb 



(i2^(cxp) 



(32) 



where dV(F O-) is the fixed-order perturbative result, dV(''='^' is the resummed cross section discussed above, and 
the truncation of the resummed cross section to the same perturbative order as d^a^^ '^ \ Here, we have 
removed the scale dependences for brevity. 

At NLO, the double-differential partonic cross section 

,(F.O.) 

(33) 



'^^'^o' ■ .(z;a„MH) = %f^)^(l-z)al°) + ^iMa«(z) + 0(a^) 



dAP dq^ ' 

receives contributions from the emission of an extra gluon jet and from processes with an initial gluon splitting into 
a qq pair, 

(34) 
(35) 
(36) 





(^) = 


Tr 

— Aqg{s,t,u) 






(^) = 


^Aqg{s,U,t) 






(^) = 


Y^Aqq{s,t,u) 





with ISC 



Aqg{s,t,u) 



s t 2u AP 



St 



(37) 



Agq{s,t,u) = -Agg{u,t,s) 



(38) 

The Mandelstam variables s, t, and u refer to the 2 — > 2 scattering process ab — > 7, + X and are related to 

the invariant mass M (or scaled squared invariant mass z = AP / s), transverse momentum qx, and rapidity y of the 
slepton pair by the well-known relations 

S = XaXbS = AP/z, (39) 



t = AP - ^JS{A/P + ql)xbey, 
u ^ AP - ^JS{AP + ql)xae-y 



(40) 
(41) 



Integration over qr requires the cancellation of soft and collinear singularities with virtual contributions in order to 
arrive at the finite single-differential partonic cross section 



da 



(F.O.) 



dM2 
,(0) 



-{z■,as,^J.R,^J.F) = aft '^(1 - z) + 



«s(M-R) ~(1) 



(42) 



where the first term ct^^ is defined in Eqs. ^ and ([8]) and the second term including the full NLO SUSY-QCD 
corrections can be found in Ref. p^ . 

The expansion of the resummed result reads 



dM2 dg2 



{N, &; as, /ifl, /if) = V fa/ ha {N + l;fJ-F) fb/h, {N + I; fip) ^ab^'' [N, b-.a^, ^ir, ^ip), 



(43) 



where a'fi^'^'^ is obtained by perturbatively expanding the resummed component 




, {N, lnx;Mfl,/iF) 



n 



(44) 
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The perturbative coefRcients E^"^ are polynomials of degree 2n in Inx, and Ti^") embodies the constant part of the 
resummed cross section in the limits 6 — > c» and N — > oo. In particular, the first-order coefficient E^^-* is given by 

SiLcc- In X) = ^a^c. X + ^^AN) In X, (45) 

with 

= -2 A(i)(5,,<5e, and sii'i)^,(iV) = -2 (s(i'<5,a(5,-, + <5,,7SI(^) + 5-cbl%{N)) . (46) 



IV. NUMERICAL RESULTS 



We now present numerical results for the production of a right-handed selectron pair at the LHC for a centre- 
of-mass energy of V5 = 14 TeV. For the masses and widths of the electroweak gauge bosons and the mass of 
the top quark, we use the values mz = 91.1876 GeV, mw = 80.403 GeV, Tz = 2.4952 GeV, Tw = 2.141 GeV, 
and mt = 174.2 GeV The electromagnetic fine structure constant a — \/^G prn^ s\v? 9w / is calculated in 
the improved Born approximation using the world average value Gp = 1.16637 • 10~^ GeV~^ for Fermi's coupling 



constant, and 



w 



m 



w 



jm 



We choose the mSUGRA benchmark point BFHK B [31|, which gives after the renormalization group evolution of 
the SUSY-breaking parameters [13] a light e/j of mass mg„ = 186 GeV and rather heavy squarks with masses around 
800-850 GeV. The top-squark mass eigenstate t\ is slightly lighter, but does not contribute to the virtual squark loops 
due to the negligible top-quark density in the proton. For the LO (NLO and NLL) predictions, we use the LO 2001 
[3^ (NLO 2004 [sl]) MRST sets of parton distribution functions. For the NLO and NLL predictions, is evaluated 
with the corresponding value of A^^^^ = 255 MeV at two- loop accuracy. We allow the unphysical scales [ip and \xr 
to vary between M/2 and 2M to estimate the perturbative uncertainty. 

In Fig. [U we present the transverse-momentum spectrum of the selectron pair, obtained after integrating the 
equations of Sec. [TTl and Sec. IIIII over M^, from the e^elj. production threshold up to the hadronic centre-of-mass 
energy. We plot the fixed order result at order ag (dashed line), the expansion of the resummed formula at the 
same perturbative order (dotted line), the total NLL+LO matched result (solid line), and the uncertainty bands 
from the scale variation. The fixed order result diverges as expected as tends to zero. The asymptotic expansion 
of the resummation formula is in good agreement with the 0{as) result in this kinematical region, since the cross 
section is dominated by the large logarithms that we are resumming. For intermediate values of g-r, we can see that 
the agreement between the expansion and the perturbative result is slightly worse. This effect was not present in 
QT-resummation for such gj^- values [l3| and is thus related to the threshold-enhanced contributions important in the 
large-M region. This can also be seen in Fig. [2l where we directly compare the jointly- and gr-niatched results, the 
latter having been obtained with the qr-resummation formalism of Ref. [25| . The two approaches lead to a similar 
behaviour in the small-gT region, but the jointly-resummed cross section is about 5%-10% lower than the gr-resummed 
cross section for transverse momenta in the range 50 GeV < qr < 100 GeV. However, the effect of the resummation 
is clearly visible in both cases, the resummation-improved result being even 40% higher than the fixed-order result 
at qT = 80 GeV. In Fig. [2 we also estimate the theoretical uncertainties through an independent variation of the 
factorization and renormalization scales between M/2 and 2M and show that the use of resummation leads to a clear 
improvement with respect to the fixed-order calculation. In the small and intermediate gT-regions the scale variation 
amounts to 10% for the fixed-order result, while it is always less than 5% for the matched result. 

The gT-distribution is affected by non-perturbative effects in the small gr-region coming, for instance, from partons 
with a non-zero intrinsic transverse-momentum inside the hadron and from unresolved gluons with very small trans- 
verse momentum. Global fits of experimental Drell-Yan data allow for different parameterizations of these effects, 
which can be consistently included in the resummation formula of Eq. (j23p through a non-perturbative form factor 
. We include in our analysis three different parameterizations of this factor [H, [s^, [s^l , 



Fab 'ib,M,xi,X2) = exp 

Fj(^^--)(6,M,x„.,) = exp 
^^r,™(6,M,xi,:.,) = exp 



-b^ ( 91+92 In ^'"^^^^'^ ]-bgig3 ln(100 a;ia;2) 



-b'^ ( 91 + 92 In ^"'''^^^ +gig3ln(100xia;2~ 



b^ \ ai + a2 In 



2 
M 



3.2GeV 



03 ln(100 xi X2] 



(47) 
(48) 
(49) 



The most recent values for the free parameters in these functions can be found in Refs. [3a, l33|- We show in the 
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FIG. 1: Transverse-momentum distribution for the process pp — > c_h e|j at the LHC. NLL+LO matched (full), fixed order 
(dotted) and asymptotically expanded results are shown, together with three different parameterizations of non-perturbative 
effects (insert). 
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FIG. 2: Transverse-momentum distribution of selectron pairs at the LHC in the framework of ioint (full) and qt (dotted) 
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upper-right part of Fig. [T] the quantity 



which gives thus an estimate of the contributions from the different NP parameterizations (LY-G, BLNY and KN). 
They are under good control, since they are always less than 5% for qx > 5 GeV and thus considerably smaller than 
the resummation effects. 

The invariant-mass distribution M^da/dM for e/j-pair production at the LHC is obtained after integrating the 
equations of Sec. |TT] and Sec. IIIII over and is shown in Fig. [31 The differential cross section da/dM has been 
multiplied by a factor in order to remove the leading mass dependence of propagator and phase space factors. We 
can see the P-wave behaviour relative to the pair production of scalar particles, since the invariant-mass distribution 
rises above the threshold at y/s — 2mg^ with the third power of the slepton velocity and peaks at about 200 GeV 
above threshold (both for M^da/dM and the not shown da/dM differential distribution), before falling off steeply 
due to the s-channel propagator and the decreasing parton luminosity. In the large-M region, the resummed cross 
section is 30% higher than the leading order cross section, but this represents only a 3% increase with respect to 
the NLO SUSY-QCD result. In the small-M region, much further then from the hadronic threshold, resummation 
effects are rather limited, inducing a modification of the NLO results smaller than 1%. The shaded, horizontally, 
and vertically hashed bands in Fig. [3] represent the theoretical uncertainties for the LO, NLO SUSY-QCD, and the 
jointly-matched predictions. At LO the dependence comes only from the factorization scale and increases with the 
momentum- fraction x of the partons in the proton (i.e. with M), being thus larger in the right part of the figure. This 
dependence is largely reduced at NLO due to the factorization of initial-state singularities in the PDFs. Including 
the dependence due to the renormalization scale in the coupling as{fiR), the total variation is about 7%-ll%. After 
resummation, the total scale uncertainty is finally reduced to only 7%-8% for the matched result, the reduction being 
of course more important in the large-M region, where the resummation effects are more important. 

In Fig. m we show the cross section correction factors 

_ daVdM 

^ ~daLO/dM ^^^^ 

as a function of the invariant-mass M. i labels the corrections induced by NLO QCD, NLO SUSY-QCD, joint- and 
threshold-resummation (as obtained in these two last calculations being matched with the NLO SUSY-QCD 

result. At small invariant mass M, the resummation is less important, since we are quite far from the hadronic 
threshold, as shown in the left part of the plot. At larger M, the logarithms become important and lead to a larger 
increase of the resummed if-factors over the fixed-order one. We also show the difference between threshold and joint 
resummations, which is only about one or two percents. This small difference is due to the choice of the Sudakov 
form factor Q and of the 7i-function, which correctly reproduce transverse-momentum resummation in the limit of 
b oo, N being fixed, but which present some differences in the pure threshold limit 6^0 and N oo, as it was 
the case for joint resummation for Higgs and electroweak boson production [2^, |2^. However, this effect is under 
good control, since it is much smaller than the theoretical scale uncertainty of about 7%. 

V. CONCLUSIONS 

With this work we complete our programme of performing precision calculations for slepton pair production at 



hadron colliders. Together with the previous papers on transverse-momentum 13[ and threshold [18i] resummation, 
soft-gluon resummation effects are now consistently included in predictions for various distributions exploiting the qt, 
threshold, and joint resummation formalisms. We found that the effects obtained from resumming the enhanced soft 
contributions are important at hadron colliders, even far from the critical kinematical regions where the resummation 
procedure is fully justified. The numerical results show a considerable reduction of the scale uncertainty with respect 
to fixed order results and also a negligible dependence on non-perturbative effects, introduced through different 
Gaussian-like smearings of the Sudakov form factors. These features lead to an increased stability of the perturbative 
results and thus to a possible improvement of the slepton pair (slepton-sneutrino) search strategies at the LHC. 
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FIG. 3: Invariant-mass distribution M^da/dM of ejj pairs at the LHC. We show the total NLL+NLO jointly matched (full), 
as well as the fixed-order NLO SUSY-QCD (dashed) and LO QCD (dotted) results, with the corresponding scale uncertainties 
(vertically hashed, horizontally hashed, and shaded bands). 
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APPENDIX A: SFERMION MIXING 



The soft SUSY-breaking terms ^/ of the trihnear Higgs-sfermion-sfermion interaction and the ofF-diagonal Higgs 
mass parameter fi in the MSSM Lagrangian induce mixings of the left- and right-handed sfermion eigenstates fL,R of 
the electroweak interaction into mass eigenstates /i.2- The sfermion mass matrix is given by [§| 



'"■LL ^ '"-/ 

TO/ rriLR 



with 



mLR 



4 

,2 



3/ sin^ 9w)m\; cos 2/3, 



m'p, + e/ sin^ 6w "m^z cos 2/3, 



At 



cot/3 for up — type sfermions. 
^ ' tan/3 for down — type sfermions. 



It is diagonahzed by a unitary matrix , M.^ S^^ — diag {rn\, m\), and has the squared mass eigenvalues 



I'LL 



'■RR 



{ml 



"RRl 



'"1,2 — '"/ ■ 2 

For real values of to^^, the sfermion mixing angle Or, < Or < 7r/2, in 



4mj |TOL_Rpj. 



cos 6 ^ sin 6 ^ 
sinOj: cos 6*^ 



with 



Jr 



can be obtained from 



tan 26 j = 



2 m/ rriLR 



to 



LL 



"-RR 



(Al) 



(A2) 
(A3) 

(A4) 



(A5) 



(A6) 



(A7) 



If tulr is complex, one may first choose a suitable phase rotation = e^'^fu to make the mass matrix real and then 
diagonalize it for /l and /jj. tan/3 = Vu/vd is the (real) ratio of the vacuum expectation values of the two Higgs 
fields, which couple to the up-type and down-type (s)fermions. The soft SUSY-breaking mass terms for left- and 
right-handed sfermions are mp and mp, respectively. 
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